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We study the ground state properties of the geometrically frustrated Hubbard model on the 
anisotropic checkerboard lattice with nearest-neighbor hopping t and next nearest-neighbor hop- 
ping t' . By using the path-integral renormalization group method, we study the phase diagram 
in the parameter space of the Hubbard interaction U and the frustration-control parameter t'/t. 
Close examinations of the effective hopping, the double occupancy, the momentum distribution and 
the spin/charge correlation functions allow us to determine the phase diagram at zero temperature, 
where the plaquette-singlet insulator emerges besides the antiferromagnetic insulator and the para- 
magnetic metal. Spin-liquid insulating states without any kind of symmetry breaking cannot be 
found in our frustrated model. 
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Strongly correlated electron systems with frustrated 
lattice structures have attracted much interest recently. 
Typical examples are a spinel compound L1V2O4 1 and 
a pyrochlore compound TI2R.U2O7 2 ' 3 , where the heavy 
fcrmion behavior and the Mott transition without mag- 
netic ordering were observed. In these compounds, elec- 
tron correlations on the frustrated lattice should play a 
vital role in yielding a variety of intriguing properties at 
low temperatures. Stimulated by the above experimental 
findings, the Hubbard model on the frustrated pyrochlore 
lattice and its two-dimensional (2D) analog, called the 
checkerboard lattice, have been studied intensively. 4-30 

In the previous paper, 20 we studied the zero- 
temperature properties of the Hubbard model on the 
isotropic checkerboard lattice at half filling, and found 
that the system undergoes a first-order phase transition 
to the plaquette-singlet insulating (PSI) phase at a finite 
Hubbard interaction. Since the analysis was focused only 
on the fully frustrated model, it is desirable to compare 
it with less frustrated models, in order to clarify how 
the frustration affects the nature of the metal-insulator 
transition. In this connection, we recall that the checker- 
board lattice is continuously connected to the ordinary 
square lattice with electron hopping t by reducing the 
amplitude of electron hopping t' along diagonal bonds 
(see Fig. 1). For the square lattice without frustration, 
it is known that the introduction of infinitesimal Hub- 
bard repulsion induces the metal-insulator transition to 
the antiferromagnetic insulating (AFI) phase. Therefore, 
the anisotropic checkerboard lattice, where the amplitude 
of diagonal hopping t' is modulated to interpolate the 
above two limiting cases, enables us to clarify the role of 
frustration in the checkerboard lattice Hubbard model. 

Motivated by this, we here investigate the Hubbard 
model on the anisotropic checkerboard lattice (Fig. [1] 
(a)) at half filling. So far, theoretical investigations in 
this direction have been put forward only in the strong 
coupling limit, where the half-filled Hubbard model can 
be mapped to the Heisenberg model with the exchange 
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FIG. 1: (color online), (a) Anisotropic checkerboard lattice. 
The solid (dotted) lines correspond to the transfer integral 
t (t 1 ). Doubly degenerate plaquette valence-bond ordering 
pattern is shown schematically. Within each unit cell the sub- 
lattice sites are denoted by 1,2. (b) The free electron band 
structure e a {k) — /io along symmetry lines in the Brillouin 
zone (B.z.) for different choices of t'/t, where ^0 is the chem- 
ical potential at half filling. 



coupling J = 4t 2 /U (</' = 4i /2 /[/). Intensive studies on 
the spin 1 /2 Heisenberg model on the anisotropic checker- 
board lattice 31-34 concluded that for large (small) J'/ J, 
the plaquette valence-bond crystal (AF Neel) phase is 
stabilized, and the corresponding phase transition is of 
first order. In this paper we aim at establishing the 
ground-state phase diagram of the anisotropic checker- 
board lattice Hubbard model in the wide parameter 
region from the weak to strong frustration/correlation 
limit. To this end, we investigate the quantum phase 
transitions at zero temperature by means of the path- 
integral renormalization group (PIRG) method devel- 
oped by Imada group, 35 " 38 which is particularly efficient 
to study electron correlations under strong frustration. 
We compute the effective hopping, the double occupancy, 
the momentum distribution and the spin/charge correla- 
tion functions. The phase diagram thus determined has 
the plaquette-singlet insulator besides the antiferromag- 
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netic insulator and the paramagnetic metal. We do not 
find spin-liquid insulating states without any kind of sym- 
metry breaking in our phase diagram. 

The paper is organized as follows. In §2, we introduce 
the model Hamiltonian and briefly explain the PIRG 
method. We discuss the quantum phase transitions to 
obtain the ground-state phase diagram of the anisotropic 
checkerboard Hubbard model in §3. A brief summary is 
given in §4. 

II. MODEL AND METHOD 

We consider the standard single-band Hubbard model 
on the anisotropic checkerboard lattice, 

^ = — tijmm ^ima^jm' a + ^ ^ ] ^im^imj i (1) 

2,^,(7, (m,m ) ? ' m 

where Ci ma (ct ) is an annihilation (creation) operator 
of an electron in the ?-th unit cell with spin a and sublat- 
tice index m (=1,2), and h lma = c\ ma c.i rncr . U is the Hub- 
bard repulsion and ty mm ' (= t, t') is the transfer integral, 
which is schematically shown in Fig. [l](a). Here, by tun- 
ing the ratio of t'/t systematically, we study the ground 
state properties to clarify how the geometrical frustration 
affects quantum phase transitions in the system. The ki- 
netic term of the Hamiltonian Ji k is diagonalized at each 
k. Then we obtain H k = J2k, a ,a £ a(k)al a(7 a kaa with two 
eigenvalues, 

£i,2(fc) = t'(cos k x + cos k y ) (2) 

±\jt' 2 ( C Osk x - COS k y ) 2 + 16i 2 (cOS 2 y)(cOS 2 5), 

where + and — signs correspond to the bands £i(k) and 
£2(fc). We assume t, t' > 0, hereafter. The dispersion 
relations thus obtained are shown in Fig. [T] (b) along the 
symmetry lines in the Brillouin zone (B.z.). 

We start by mentioning some characteristics in the two 
limiting cases of t'/t = and t'/t = 1. In the case of 
square lattice (t'/t = 0), it is known that the AFI state 
is stabilized for any finite U > at zero temperature, due 
to the perfect nesting. In our reduced B.z. scheme, the 
nesting is expressed by the hybridization between upper 
and lower bands with the same momentum along the 
boundary of B.z. as shown in Fig. [1] (b). On the other 
hand, for the isotropic checkerboard lattice (t'/t = 1), 
the upper band is completely flat over the whole B.z. 
while the lower band is dispersive at U — 0. As shown 
in Refs. 14,15 , the perturbative calculation in U at half 
filling in the isotropic case suffers from divergence at third 
and higher orders, because the lower band is completely 
filled and the Fermi level just touches the flat band. This 
unusual situation makes the theoretical treatment of the 
model difficult. 

For highly frustrated lattice systems, it is known that 
powerful quantum Monte Carlo method suffers from the 



minus sign problem. Also, the exact diagonalization cal- 
culation cannot deal with large enough lattice sizes to fig- 
ure out the role of frustration in our system. To treat the 
strong correlation and frustration effects, we here make 
use of the PIRG method 35-37 , where we further employ 
a quantum number projection (QP) operator to the total 
spin-singlet state. 37 The PIRG+QP method is particu- 
larly efficient to study electron correlations under strong 
frustration. The algorithm is very simple. We start from 
an unrestricted Hartree-Fock solution 39 and reach the 
ground state by taking into account quantum fluctua- 
tions in a systematic fashion. We increase the dimension 
of the truncated Hilbert space in a nonorthogonal Slater 
basis numerically optimized by the path-integral opera- 
tion. An energy variance extrapolation 35-38,40 is very 
efficient to reach the true ground state of finite size sys- 
tems. We take the number of Slater basis functions up to 
500 and apply an improved iteration scheme proposed in 
our previous paper 20 . In the present study, we carry out 
the calculation for the N = 32 lattice system with pe- 
riodic boundary conditions. We have already confirmed 
that the N = 32 lattice system is large enough to inves- 
tigate the quantum phase transitions in the thermody- 
namic limit in the isotropic case 20 . 

III. RESULTS 

To study the quantum phase transitions at zero tem- 
perature, we first consider the virtual displacement of E g 
with respect to t' and U, 

tE^.m-^u+^a, (3) 

and evaluate each coefficient. The former coefficient de- 
scribes the averaged hopping amplitude along diagonal 
bonds, while the latter the double occupancy of electrons 
at each site. Both quantities provide important informa- 
tion about the quantum phase transition. 

To discuss the spin properties, we also calculate the 
site-dependent spin correlation function C t ( t <) defined by 

N N t 

c * = iiEE( s <'^). ( 4 ) 

j=l T t = l 

where r t labels neighboring sites connected by the trans- 
fer integral t and N t is the number of them (Nt = 4, 
Nf = 2). Cf is obtained by replacing r t with t# in the 
formula (4). For our anisotropic checkerboard lattice, the 
plaquette spin singlet state is one of the most probable 
candidates for the ground state, so that we examine the 
following plaquette correlation function Pa(b) defined by 

Pa = (Q\), 

9 N/2 r*\ 

Qa = 4 E ( } 

i=l 
pattern A 
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where pi [= (S ai + S 7i )-(Sp i + Ss t )] is the ith plaquette- 
singlet operator. The corresponding configuration pat- 
tern of plaquettes and their signs (—1)* are schematically 
shown in Fig. Q] (a), where filled (open) circles represent 
positive (negative) signs for Pa- The pattern for the cor- 
relation Pb is given by totally shifting the circles to the 
square with crossing in Fig. [T] (a). Note here that in 
the autocorrelation terms X)j(pf) m Eq. (5) equally con- 
tribute to Pa and Pb in the AF ordered phase, so that 
the plaquette correlation functions Pa and Pb could have 
finite values even in this phase. Therefore, we also ex- 
amine the value of Pa — Pb in order to distinguish the 
plaquette singlet phase and the AF phase. 
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FIG. 2: (color online). (a) The averaged hopping, 
— (dE g /dt') u , (b) the spin correlation CW), and (c) the pla- 
quette correlation Pa(b) as a function of t'/t on the TV = 32 
lattice system at half filling for different choices of U/t. The 
inset of (c) shows Pa — Pb ■ 



A. Phase transitions under frustration control 

In the following, we present the computed results in 
two ways; quantum phase transitions are discussed un- 
der control of (i) frustration and (ii) electron correlations. 
In this subsection, we first address the frustration con- 
trol, and show the above physical quantities as a func- 
tion of t'/t for different choices of U/t. The results are 
summarized in Fig. [21 Starting from the square lattice 
with t' = 0, we examine the instability of the AFI state 
which is stabilized at t' = for any finite U/t > 0. In 
the AFI phase, two electron spins sitting on the nearest 
neighbor (next nearest neighbor) sites have AF (ferro- 
magnetic) correlations, as indeed confirmed in Fig. [2Kb). 
Correspondingly, the absolute value of the averaged hop- 
ping on the t' bond, (dE g /dt') u , is strongly suppressed 
in the AFI phase (Fig. [5] (a)). These results indicate 
the presence of AF order in the small t'/t region at least 
for the Hubbard repulsion, U > 4. We also note that 
the double occupancy (dE g /dt') u is almost unchanged 
in the AFI phase even if t'/t is altered (though not shown 
in the figure). When the effect of frustration is further 
enhanced via increase of t'/t, the AFI state becomes un- 
stable, triggering a phase transition to another insulating 
phase. The transition point depends on the strength of 
the Hubbard U. We indeed observe the abrupt jump in 
(dE g /dt') u at t'/t ~ 0.4, 0.6, and 0.8 for U/t = 4, 6, 
and 8, respectively (Fig. [2(a)). Similar discontinuities in 
the two types of correlation functions are found at the 
same transition points in Figs. [2| (b) and (c). Therefore 
we conclude that the phase transition is of first order. 

In order to see the nature of the insulating phase re- 
alized at larger t'/t in detail, we look at the spin cor- 
relation C t (f) and the plaquette correlation Pa(b)- In 
the case of U/t = 4 and 6, we find that the first-order 
transition occurs between AFI and PM, as clearly seen in 
C* t (t') shown in Fig. [21 (b). In these cases, once the tran- 
sition occurs Pa and Pb are both reduced and the value 
of Pa — Pb is almost zero, implying that the plaquette- 
singlet state is not formed. We also note that in the PM 
phase — (dE g /dt') u has the value about 0.3, which is a 
little bit smaller than 1 /3 expected for the noninteracting 
case (Fig. (21(a)). 

In contrast to the above two cases, quite different be- 
havior emerges in the correlation function in the case of 
U/t = 8: not only the magnitude of Pa(b) itself but also 
Pa — Pb are abruptly increased at t = t' c , as seen in Fig. 
[21(c) and its inset, implying that the transition from AFI 
to PSI indeed occurs. We have confirmed that the same 
type of transition occurs at t'/t ~ 0.9 for U/t = 10. In 
the strong coupling limit with large U, we can check how 
precise our estimate of the transition point is in terms of 
the effective model. In this limit, our system is mapped 
to the anisotropic checkerboard Heisenberg model with 
two different exchange couplings J and J'. The value of 
t'Jt for the transition point then yields J'J J ~ 0.8, which 
agrees very well with J' c /J = 0.79 — 0.81 for the Heisen- 
berg model estimated by a strong-coupling expansion 33 . 
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The good agreement confirms the validity of our analysis, 
and in turn supports the existence of the first order phase 
transition between AF Neel and plaquette valence-bond 
crystal phase proposed for the Heisenberg model 32-34 . 
We note here that the above results are totally consistent 
with our previous study of the isotropic model (t'/t = 1), 
where much more systematic analyses, performed with 
finite-size scaling, give U/t = 6.75 ± 0.25 20 for the PM- 
PSI transition point 20 . 



Phase transitions under control of electron 
correlations 
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FIG. 3: (color online). (a) The double occupancy 
(dE g /dU) t ,, (b) the spin correlation CW), and (c) the pla- 
quette correlation -Pa(b) as a function of U/t on the N = 32 
lattice system at half filling for different choices of t' /t. The 
inset of (c) shows Pa — Pb- 



We next show the physical quantities in Eqs. (3)-(5) 
as a function of U/t for different choices of t'/t to discuss 
the quantum phase transitions under control of electron 
correlations. We show the [7/i-dependence of the dou- 



ble occupancy (dE g /dU) t , in Fig. [3] (a). The introduc- 
tion of the Hubbard interaction monotonically decreases 
the double occupancy, implying that the paramagnetic 
metallic state is realized in the small U region. Fur- 
ther increase in the interaction yields the discontinuity 
of (dE g /dU) t ,, in accordance with the first order Mott 
transition. We determine the transition point U c by es- 
timating the level crossing point of energies for the com- 
peting metallic and insulating states. 

In the case of t'/t = 0.5, C t and C t > both show disconti- 
nuities at U c /t ~ 5.4 and suddenly increase with opposite 
signs, suggesting the transition to the AFI phase. For 
larger U (> U c ), the absolute value of C t is smaller than 
that for Ct — —0.335 known for the square lattice Heisen- 
berg model (J' / J = 0) 41 . The difference comes from the 
presence of charge fluctuations which reduce the local 
magnetic moment. For t'/t = 0.8, the phase transitions 
take place twice at U c \/t ~ 6.6 and U c i/t ~ 7.6, and the 
PSI phase is realized for U c \ < U < U c % between the PM 
and AFI phases (see Fig. [3](c) and its inset). Therefore 
we can see three types of phase transitions under con- 
trol of correlations: a single Mott transition of PM-AFI 
(t'/t = 0.5) and PM-PSI (t'/t = I), and double quantum 
phase transitions of PM-PSI and PSI- AFI (t'/t = 0.8). 

To investigate the nature of the quantum phase transi- 
tions in detail, we further calculate the momentum distri- 
bution n a (k) and the momentum-dependent correlation 
functions in the charge [spin] sector N mm ' (q) [S mm > (q)] 
at t'/t — 0.8, which are given by 



n a (k) 



^^2( & La & kaa) for fc ^ (tT, 7t), 

^E^U^) for k = 7r )> 



(6) 



a,0 



„ N/2 

2 . . / \ 

i,3=l 



x e 



N/2 



x e 



\q-{Ri m -R,i) 



(7) 



(8) 



where hi m — hi m ^ + fiimi and Ri m represents the posi- 
tion of the i-th unit cell in the m-th sublattice. Diag- 
onalizing the 2x2 matrix, we obtain Sp(q) and Np(q) 
(f3 = max, min). 

We show the computed results in Fig. 4. In the PM 
phase, the quasi- Fermi surface exists and therefore n%(k) 
has a discontinuity at k — (tt, tt) for U < U c \ (Fig. 4 
(a)). As U increases, the PM state becomes unstable and 
then the discontinuity of ri2(fe) disappears at U > U c \ in 
the insulating phase (Fig. 4 (b)-(c)). Moreover we can 
confirm that the Hubbard interaction reduces the charge 
fluctuations so that Np (q) is uniformly reduced over the 
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FIG. 4: (a)-(c) The momentum distribu- 
tion function n a (k), where a = 1 (2) cor- 
responds to the dashed (solid) lines and 
(d)-(f) [(g)-(i)] the momentum-dependent 
charge [spin] correlation function Np(q) 
[S/s(q)], where /3 =max (min) corresponds 
to the solid (dashed) lines at t'/t = 0.8 
for different choices of U/t. Note that 
4Sp(q)=N {q) at U/t = 0. 
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whole B.z. (Fig. 4 (d)-(f)). By contrast, the spin corre- 
lations 5 max (g) are totally enhanced. Especially a peak 
structure is developed for U > U c \ at q — (0, 0) (Fig. 4 
(h)-(i)), while there is no such peak structure for U < U c i 
(Fig. 4 (g)). Also, divergent increase in S max (0, 0) is ob- 
served beyond a certain interaction U c %. Therefore the 
phase transition between PSI and AFI occurs at U = U C 2- 
As a result of the transition, the momentum distribu- 
tions ni(k) and n2(fe) have the reflection symmetry with 
respect to the n a (k) — 0.5 plane (Fig. 4 (c)). These 
properties in the AFI phase should be adiabatically con- 
nected to those in the SDW phase for the square lattice 
Hubbard model (t'/t = 0). 

Summarizing all the above results, we end up with the 
phase diagram of the Hubbard model on the anisotropic 
checkerboard lattice, as shown in Fig. [5j There are three 
distinct phases of PM, AFI and PSI. The quantum phase 
transitions among them are of first order. For small U, 
the PM phase appears, while for large U and small t', 
the AFI phase is realized in accordance with the known 
results. The PSI phase is stabilized in the strong frus- 
tration region with t' jt ~ 1 and large U. Note that 
around t' jt ~ 0.8 and U/t ~ 6, the three phases strongly 
compete with each other. In fact, the double phase tran- 
sitions, which occur around t'/t = 0.8 as U increases, 
reflect this kind of strong competition. We would like to 
stress again that in the two limiting cases with strong 
frustration, the phase boundary obtained here reproduce 



the known results fairly well: the Heisenberg limit 33 
with large U and the isotropic checkerboard limit 20 with 
t'/t = 1, both of which were studied in detail previously. 
We therefore believe that the phase diagram obtained 
in this paper is reliable although we have restricted our 
analysis to the TV = 32 lattice system. 



IV. SAMMARY 

We have studied the ground-state properties of the 
anisotropic checkerboard lattice Hubbard model by 
means of the PIRG method. By controlling the geomet- 
rical frustration via a systematic change in the transfer 
integral t' along diagonal bonds, we have dealt with the 
wide parameter region from the square lattice {t! = 0) to 
the fully- frustrated isotropic checkerboard lattice (t'/t = 
1). The ground state phase diagram thus obtained con- 
sists of three distinct quantum phases depending on the 
Hubbard interaction U and the strength of frustration 
t'/t. 

In particular, in the region with strong frustration 
(t'/t ~ 1), we have the plaquette singlet insulator 
with broken translational symmetry, which is in con- 
trast to the results for analogous two-dimensional frus- 
trated electron systems such as the anisotropic triangular 
lattice model where a quantum spin liquid phase with- 
out any symmetry breaking was proposed for the ground 
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state 42 ~ 49 . Therefore, it remains an important problem 
to figure out what is really relevant for realizing the insu- 
lating phase without symmetry breaking. It is also inter- 
esting to investigate the nature of the finite-temperature 
Mott transition of the anisotropic checkerboard lattice. 
In particular, it is worth exploring whether the reen- 
trant behavior in the temperature-driven Mott transition 
found for the anisotropic triangular lattice model 50 could 
emerge in the checkerboard lattice model. These issues 
are now under consideration. 
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FIG. 5: The ground-state phase diagram of the anisotropic 
checkerboard lattice Hubbard model. AFI, PSI, and PM rep- 
resent the antiferromagnetic insulator, the plaquette singlet 
insulator and the paramagnetic metal, respectively. Calcu- 
lated results for the N = 32 system are shown by open circles 
with error bars. The phase transition point is determined by 
the crossing point of the energies for two competing phases, 
and its error bars originate from the extrapolation procedure 
of the energies. Cross symbol indicates the phase boundary 
of the U — » oo Heisenberg limit 33 (see text). The quantum 
phase transitions are of first-order. 
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